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Abstract 

A combination of small-cluster exact-diagonalization calculations and a 
well-controlled approximative method is used to study the ground-state phase 
diagram of the spin-one-half Falicov-Kimball model extended by the spin- 
dependent on-site interaction between localized (/) and itinerant (d) electrons. 
Both the magnetic and charge ordering are analysed as functions of the spin- 
dependent on-site interaction (J) and the total number of itinerant {N^) and 
localized (Aj) electrons at selected U (the spin- independent interaction be- 
tween the / and d electrons). It is shown that the spin-dependent interaction 
(for Nf = L, where L is the number of lattice sites) stabilizes the ferromagnetic 
(F) and ferrimagnetic (FI) state, while the stability region of the antiferromag- 
netic (AF) phase is gradually reduced. The precisely opposite effect on the 
stability of F, FI and AF phases has a reduction of Nf. Moreover, the strong 
coupling between the / and d-electron subsystems is found for both Nj = L 
as well as Nf < L. 

PACS nrs.:75.10.Lp, 71.27.+a, 71.28.+d, 71.30.+h 
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1 Introduction 



In the past decade, a considerable amount of effort has been devoted to understand 
the underlying physics that leads to a charge ordering in strongly correlated electron 
systems. The motivation was clearly due to the observation of a such ordering in 
doped nickelate P and cuprate |2] materials, some of which constitute materials that 
exhibit high-temperature superconductivity. One of the simplest models suitable to 
describe charge ordered phases in interacting electron systems is the Falicov-Kimball 
model (FKM) Indeed, it was shown that already the simplest version of this 
model (the spinless FKM) exhibits an extremely rich spectrum of charge ordered 
solutions, including various types of periodic, phase-separated and striped phases 
However, the spinless version of the FKM, although non-trivial, is not able to account 
for all aspects of real experiments. For example, many experiments show that a 
charge superstructure is accompanied by a magnetic superstructure 12] • In order 
to describe both types of ordering in the unified picture Lemanski [S] proposed a 
simple model based on a generalization of the spin-one-half FKM with an anisotropic, 
spin-dependent local interaction that couples the localized and itinerant subsystems. 
The model Hamiltonian is 



= X! ^ijdtadja + fi^fiadia'dia' + J '^{ft-afi~a " filjfia)dtadia, (1) 

ijcT i(T(T' i(T<T' 

where f^, fia are the creation and annihilation operators for an electron of spin 
cr =t, i in the localized state at lattice site i and df^, di^ are the creation and 
annihilation operators of the itinerant electrons in the rf-band Wannier state at site i. 

The first term of (1) is the kinetic energy corresponding to quantum-mechanical 
hopping of the itinerant d electrons between sites i and j. These intersite hopping 
transitions are described by the matrix elements tjj, which are —t if i and j are the 
nearest neighbors and zero otherwise (in the following all parameters are measured 
in units of t). The second term represents the on-site Coulomb interaction between 
the (i-band electrons with density = N^/L = jJ2ia df^di^- and the localized / 
electrons with density nj = Nf /L = J2ia f tafias where L is the number of lattice 
sites. The third term is the above mentioned anisotropic, spin-dependent (of the 
Ising type) local interaction between the localized and itinerant electrons that reflects 
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the Hund's rule force. Moreover, it is assumed that the on-site Coulomb interaction 
between / electrons is infinite and so the double occupancy of / orbitals is forbidden. 

Thus from the major interaction terms that come into account for the interacting 
d and / electron subsystems only the Hubbard type interaction between the spin-up 
and spin-down d electrons has been omitted in the Hamiltonian (1). In his work j3] 
Lemanski presents a simple justification for the omission of this term, based on an in- 
tuitive argument: the longer time electrons occupy the same site, the more important 
becomes interaction between them. According to this rule the interaction between 
the itinerant d electrons {Udd) is smaller than the interaction between the localized 
/ electrons (f///) as well as smaller than the spin-independent interaction between 
the localized and itinerant electrons U . In this paper we specify more precisely con- 
ditions when this term can be neglected. For this reason we start our study with the 
case Udd 7^ 0. To determine the effects of Udd interaction on the ground-states of the 
conventional spin-one-half FKM ( J = 0) the exhaustive studies of the ground-state 
phase diagrams of the model (in the nj — Udd plane) are performed for several cluster 
sizes. Of course, an inclusion of the Udd term makes the Hamiltonian (1) intractable 
by methods used for the conventional sp in-one-half /spinless FKM and thus it is nec- 
essary to use other numerical methods. Here we use the Lanczos method to study 
exactly the ground states of the spin-one-half FKM generalized with Udd interaction 
between the spin-up and spin-down d electrons. 

2 Results and discussion 

2.1 The spin-1/2 FKM with the Hubbard interaction be- 
tween itinerant electrons 

The Hamiltonian of the spin-1/2 FKM with the Hubbard interaction between itin- 
erant d electrons can be written as the sum of three terms: 

H = Y. UjdtJja + ftJiadl'dia' + Udd E dt^d,^d+d,^. (2) 

ijtj itjcr' i 

Since the /-electron density operators fj^fia of each site i commute with the Hamil- 
tonian (2), the /-electron occupation number is a good quantum number, taking 
only two values, Wjo- = 0, 1 according to whether the site i is unoccupied or occupied 
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by the localized / electron of spin a. Therefore the Hamiltonian (2) can be rewritten 

as 

^ = I] hijdf^dj^ + Udd J2 df^di^df^dii, (3) 

ija i 

where hij = tij + U{wii+Wii)6ij. This Hamiltonian can be considered as a generalized 
Hubbard model. To determine the ground-state energy of the model we have used 
the Lanczos method 0- However, since the hopping amplitudes depend now on the 
/-electron distribution w = {wi, W2 ■ ■ ■ wl} (we remember that the double occupancy 
of / orbitals is forbidden Wi = + = 0, 1) the Lanczos procedure has to be 
used many times (strictly said L\/{{L — Nf)\Nf\) times) for given L and Nf. Of 
course, such a procedure demands in practice a considerable amount of CPU time, 
which imposes severe restrictions on the size of clusters that can be studied within 
the exact-diagonalization method. For this reason we were able to investigate exactly 
only the clusters up to L = 12. Fortunately, it was found that in some parameter 
regimes the ground-state characteristics of the model are practically independent of 
L and thus already such small clusters can be used satisfactorily to represent the 
behavior of macroscopic systems. In particular, we have studied the stability of the 
ground-state configuration w^{Nf) (obtained for Udd = and fixed A'^^) at finite 
values of Udd- The results of numerical calculations obtained for U = A and U = 8 
are summarized in Fig. 1 in the form of — Udd phase diagrams (the half- filled band 
case n f + Ud = 1 is considered). One can see that the ground-state configuration 
w^{Nf) found for Udd = persists as a ground state up to relatively large values of 
Udd (t^drf ~ 2, for f/ = 4 and f/J^ ~ 6, for U = 8), revealing small effects of the Udd 
term on the ground state of the model in the strong U interaction limit. Contrary to 
the strong coupling case, for small [U = 1) and intermediate {U = 2) values of the 
Coulomb interaction between the localized and itinerant electrons very strong effects 
of Udd term on the ground states of the model have been observed. In these cases the 
typical values of f/J^ are of order 0.5 but for some Nf even much smaller values were 
found. Thus we can conclude that the Hubbard type interaction between the spin-up 
and spin-down d electrons can be neglected in the strong interaction limit between 
the localized and itinerant electrons {U > 4). For this reason all next calculations 
on the spin-one-half FKM with spin-dependent Coulomb interaction J between the 
/ and d electrons have been done at [/ = 4. 
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2.2 The spin-1/2 FKM with the spin-dependent interaction 
between itinerant and locahzed electrons 

The spin-dependent interaction term Hj between the itinerant {d) and locahzed (/) 
electrons does not violate the condition [f^fiai H] = and thus the Hamiltonian (1) 
can be rewritten as 

^ = E hijd+d.^, (4) 

ijcr 

where hij = tij + {Uwi + Jwii — Jwi^)5ij. Thus for a given /-electron configuration 
w = {wi,W2 ■ ■ -Wl}, the Hamiltonian (4) is the second-quantized version of the 
single-particle Hamiltonian h{w), so the investigation of the model (4) is reduced to 
the investigation of the spectrum of h for different configurations of / electrons. This 
can be performed exactly, over the full set of /-electron distributions (including their 
spins), or approximatively, over the reduced set of /-electron configurations. The 
second way has been used in the original work by Lemanski j^j. He studied the two- 
dimensional version of the model using the method of restricted phase diagrams (all 
possible configurations of the localized / electrons for which the number of sites per 
unit cell is less or equal to 4 are considered) and presented some preliminary results 
concerning the charge and magnetic order in the ground-state of this model. For 
example, he detected various phases with complex charge and magnetic structures 
that form consecutive stages of transformation of F to AF phase with an increase 
of the band filling. In the present work we study the one (D=l) and two (D=2) 
dimensional analogue of the model. To determine the ground states of the model we 
use the method of small cluster-exact diagonalization calculations in a combination 
with a well-controlled numerical method [7j. 

Since the d electrons do not interact among themselves, the exact numerical 
calculations on finite clusters precede directly in the following steps: (i) Having U, J 
and w = {wi, W2 ■ ■ ■ wl} fixed, find all eigenvalues \k of h{w). (ii) For a given Nf = 
J^iWi and Nd determine the ground-state energy E{w) = Z]£i -^fc of a particular 
/-electron configuration w by filling in the lowest one-electron levels (the spin 
degeneracy must be taken into account), (iii) Find the for which E{w) has a 
minimum. Repeating this procedure for different values of Nf, Nd, U and J, one 
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can immediately study the ground-state phase- diagrams of the model in different 
parameter spaces. 

To reveal an influence of the anisotropic, spin-dependent interaction between the 
localized and itinerant electrons on the ground states of the model we have started 
the study with the case A'^^ = L (D=l). In this case each lattice site is occupied 
by one (up or down) / electron (the double occupancy is forbidden) and thus only 
distributions over different spin configurations should be examined. Although the 
total number of configurations increases very rapidly with the cluster size L (as 2^), 
relatively large lattices can be reached by this method {L ~ 32) if all symmetries 
of the Hamiltonian are considered. In Fig. 2 we summarize numerical results ob- 
tained by small-cluster exact diagonalization calculations on the largest cluster that 
we were able to consider exactly (L = 32) in the form of — J phase diagram. To 
avoid an ambiguity in determination of FI and AF phases we examined the ground 
states only for even Nd- In the J direction the calculations have been done with 
step AJ = 0.05. Various phases that enter into the phase diagram are classified 
according to Sj = J^i'^ii — "^ii and = Ndi — Nai'. the ferromagnetic phase is 
characterized by 15*11 = A'^^, {S^l = N^, the ferrimagnetic phases are characterized 
by < 15*11 < A^/,0 < |5J| < A^^ and the antiferromagnetic phases are character- 
ized by 15*11 = 0, |5J| = 0. Comparing numerical results obtained for |5*|| and 15*11 
one can find a nice correspondence between the magnetic phase diagrams of local- 
ized (/) and itinerant (d) subsystems. Indeed, with the exception of several isolated 
points at J = 0.05, the corresponding F, FI and AF phases perfectly coincide over 
the remaining part of diagrams showing on the strong coupling between the mag- 
netic subsystems of localized and itinerant electrons for nonzero values of J. It is 
interesting that already very small changes of the spin-dependent interaction can 
produce so important cooperative changes. This confirms the supposition that the 
spin-dependent interaction between the localized and itinerant electrons could play 
an important role in description of ground state properties of the generalized FKM. 
In general, the spin-dependent interaction J stabilizes the F and FI phases while the 
AF phase is gradually suppressed with increasing J. Moreover, we have observed 
that within the AF phase (with the exception of cases A^^^ = 14, 26) the ground states 
(for given A^^) do not change when J increases, while within the FI phase very strong 
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effects of J on ground states were found (see Fig. 2). For example, the transition 
from the AF to F phase at Nj = 12 reahzes through the following sequence of FI 
phases: [T2i6]4 -^T2ii2T2i3i2T2^T26i2T2i2-^T3oi2, where the lower index denotes 
the number of consecutive sites occupied by up or down spin / electrons, or the 
number of repetitions of the block [...]. In Fig. 2 we present also the complete set of 
ground-state configurations (obtained on the above specified set of Nd and J values) 
from the AF region. Between them one can find different types of periodic and non- 
periodic configurations, but the most interesting examples represent configurations 
formed by antiparallel F domains, that illustrate convincingly the cooperative effects 
of spin-dependent interaction J between the localized and itinerant electrons. 

Let us turn our attention to the case Nf ^ L. From the numerical point of 
view this case is considerably exacting, since now we have to minimize the ground 
state energy not only over all different spin configurations but also over all different 
/-electron distributions. This takes a considerable amount of CPU time and for 
this reason we were able to investigate exactly only the clusters up to L = 24 for 
Nf 7^ L. In Fig. 3 we present the one-dimensional skeleton phase diagram of the 
generalized FKM in the Nf—N^ plane obtained by small-cluster exact diagonalization 
calculations for L = 24, U = 4 and J = 0.5. Again, the stability regions of AF, F 
and FI phases are marked by (■), (+) and (o). Here we displayed the numerical 
results only for the localized subsystem since the analysis of \S^\ and 15*11 showed 
that the F, FI, and AF phases corresponding to localized and itinerant subsystems 
coincide also for Nf < L (similarly as for Nf = L, only a few exceptions have 
been observed in isolated points that result probably from the finite-size effects). 
The most striking feature of the A'^^ — N^ phase diagram is that with decreasing 
Nf the AF phase is stabilized, while the F and FI phases are suppressed. It is 
interesting that this effect is strongly asymmetric and a disappearance of F and 
FI phases realizes in two different ways. For small rf-electron concentrations the F 
phase disappears practically immediately outside the point Nf = L. The FI phase 
survives along the main diagonal in the narrow band and fully disappears near the 
point Nf = L/2,N(i = L. In the opposite limit (large d-electron concentrations) the 
F phase persists for a wide interval of Nf < L and with decreasing Nf disappears 
gradually. The same behaviour exhibits also the FI phase. The F phase appears also 
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for Nf — > 0. However, in this limit strong finite-size effects have been observed, and 
thus it is questionable if this phase persists really in the thermodynamic limit. 

Although the skeleton phase diagram of the generalized FKM is rather simple, 
the spectrum of magnetic solutions that yields the model for the AF and FI phases 
is very rich. Indeed, for L = 24 and J = 0.5 we have found 140 different AF phases 
and 20 different FI phases that enters into the Nf — Nd phase diagram. Of course, 
it is not possible to present here all ground-state configurations, but let us show at 
least several main configurations types, with the largest stability regions. From the 
AF phases the largest stability region (denoted by I) corresponds to configurations 
of the type t„J.„ 0L-2n- The second largest region (denoted by II) corresponds to AF 
configurations of the type [| 0„ | 0„]feOL_2fc(r!,+i)- Typical examples of the AF ground 
states from the central region of the phase diagram represent periodic configurations 
of the type tn Om in Om (below the main diagonal) and configurations of the type 
T2 [iT]fci [2 Om.[l Op t Op]fc2, or t2 [itlfci i2 0„[i Op t 0p-i]fc2, above the main diagonal. 
Between these configurations and the F region the ground states are the segregated 
configurations of the type U [iT]fc l2 0™,, or [ti]fci U [iT]fc2 U [Tilfe^Om, or their 
modifications (the region denoted by III). In the FI region the typical examples of 
ground states represent configurations of the type tn [0 i t]fcO J. 0. 

These results show that the spectrum of magnetic and charge solutions that 
yields the spin-one-half FKM model generalized with the spin-dependent interaction 
between the localized and itinerant electrons is indeed very rich. Of course, one can 
ask if these phases persist also on larger clusters. Unfortunately, lattices larger than 
L = 32 (for Nf = L), or L = 24 (for Nf < L) are beyond the reach of present 
day computers within the exact diagonalization technique. Therefore, to resolve this 
problem one has to use other numerical methods. Very promising seems to be the 
well-controlled numerical method that we have elaborated recently to study ground 
states of the spinless FKM [7]. This method is described in detail in our previous 
papers [71 |H] and thus we summarize here only the main steps of the algorithm that 
is a simple modification of the exact-diagonalization algorithm described above: (i) 
Chose a trial configuration w = {wi,W2, ■ ■ ■ ,wl}. (ii) Having w, U and J fixed, 
find all eigenvalues of h{w). (iii) For a given A'^^ = J^i'^i and N^ determine the 
ground-state energy E{w) = J2k=i -^k of a particular /-electron configuration w by 
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filling in the lowest A'^^ one-electron levels, (iv) Generate a new configuration w' by 
moving a randomly chosen electron to a new position which is chosen also at random, 
(v) Calculate the ground-state energy E{w'). If E{w') < E{w) the new configuration 
is accepted, otherwise w' is rejected. Then the steps (ii)-(v) are repeated until the 
convergence (for given parameters of the model) is reached. Of course, one can move 
instead of one electron (in step (iv)) simultaneously two or more electrons, thereby 
the convergence of the method is improved. Indeed, tests that we have performed 
for a wide range of the model parameters showed that the latter implementation 
of the method, in which 1 < p < Pmax electrons {p should be chosen at random) 
are moved to new positions overcomes better the local minima of the ground state 
energy. In this paper we perform calculations with Pmax = Nf. The main advantage 
of this implementation is that in any iteration step the system has a chance to 
lower its energy (even if it is in a local minimum), thereby the problem of local 
minima is strongly reduced (in principle, the method becomes exact if the number of 
iteration steps goes to infinity). On the other hand a disadvantage of this selection 
is that the method converges slower than for Pmax = 2 and Pmax = 3. To speed up 
the convergence of the method (for Pmax = Nf) and still to hold its advantage we 
generate instead the random number p (in step (iv)) the pseudo-random number p 
that probability of choosing decreases (according to the power law) with increasing 
p. Such a modification improves considerably the convergence of the method. 

To test the convergence of the method we have first calculated the ground-state 
configurations of the model in the Nf — A^^^ plane on the cluster of L = 24 sites. 
Comparing these results with exact ones (discussed above) we have found that the 
method is able to reproduce the exact ground states after relatively small number of 
iterations (typically 5000-10000 iterations per site). Then we have used the method 
to study the A'^^ — phase diagram of the model on larger clusters consisting of 
L = 36 and L = 48 sites. Our numerical computations showed that all main results 
obtained on small clusters hold also on larger clusters. Again we have observed 
strong coupling between two magnetic subsystems and a coincidence of corresponding 
magnetic phases, that stability regions are practically unchanged with increasing L 
(see Fig. 4). Moreover, we have observed that the main configurations types found 
for L = 24 persist also on large clusters and thus we suppose that the real magnetic 
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phase diagram of the model will be very close to ones presented in Fig. 3 and Fig. 4, 
Of course, it is possible that some ground states that are uniform on the finite lattices 
could be degenerated in the thermodynamic limit. 

The same calculations we have performed also in two dimensions. The two- 
dimensional results are of particular importance since they could shed light on the 
mechanism of two-dimensional charge and magnetic ordering in doped nickelate 
and cuprate |2] materials. Here we concern our attention on a description of basic 
types of charge and magnetic ordering that exhibits the spin-one-half FKM with spin- 
dependent interaction between d and / electrons in two dimensions. To minimize 
the finite-size effects the numerical calculations have been done on three different 
clusters of 4 x 4, 6 x 6 and 8x8 sites. On the 4x4 cluster the calculations were 
performed by exact-diagonalization method and on larger clusters the approximative 
method described above was used. 

Similarly as in the one dimension we start our two-dimensional studies with the 
case Nf = L. The magnetic phase diagrams of the / and d electron subsystems 
calculated for Nf = L [L = 36) are shown in Fig. 5. Comparing these phase 
diagrams with their one-dimensional counterparts one can find obvious similarities. 
In both cases the basic structure of the phase diagram is formed by three large F, 
FI and AF domains that are accompanied by secondary phases. However, while in 
the one dimension the secondary phases are stable only in isolated points at very 
small values of J, in two dimensions these secondary phases persist also for large J. 
Calculations that we have performed on different clusters (4 x 4, 6 x 6, and 8x8) 
showed that the secondary structure depends very strongly on the cluster-size (with 
increasing L it is gradually suppressed) and it is not excluded that it fully disappears 
in the thermodynamic limit {L — * oo). 

The typical examples of the ground-state configurations (that represent the most 
frequently appearing types of the ground states in the — J phase diagram) are 
displayed in Fig. 6. Again one can see that the spectrum of magnetic solutions that 
yields the FKM extended by spin-dependent interaction is very rich. In addition to 
the F phase (that the stability region shifts to higher ci-electron concentrantions when 
J increases) there are various types of AF and FI structures like the antiparallel F 
domains (2-3), the axial magnetic stripes (4-7), the diagonal magnetic stripes (8-11) 
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and the perturbed diagonal magnetic stripes (12-15). This again demonstrates strong 
effects of the spin-dependent interaction on the formation of magnetic superstructures 
in the extended FKM and its importance for a correct description of correlated 
electron systems. 

From the experimental point of view the most interesting case is, however, the 
case Nf < L that could model the real situation in doped nickelate and cuprate 
systems 0I21- To describe possible charge and magnetic orderings for Nf < L we 
have performed an exhaustive studies of the model on the 6x6 and 8x8 clusters. 
Numerical calculations have been done over the full set of Nf and A'^^ values and 
they revealed a rich spectrum of coexisting charge and magnetic superstructures. 
Some typical examples of these superstructures are presented in Fig. 7. Between 
them one can find various types of phase segregated (e.g., 1), phase separated (e.g., 
16) and n-molecular (e.g., 3) configurations with F,FI and AF ground states as well 
as various types of axial (e.g., 9) and diagonal (e.g., 5) magnetic /charge stripes. In 
generally, we have observed that the system has tendency towards phase segregation 
for small and large ci-electron concentrations, while near the na = 1 point the system 
prefers to form the various types of axial and diagonal stripes. In addition, similarly 
as in D = 1 a strong reduction of F and FI phases with decreasing Nf is observed 
also in D = 2. We have found that the same tendencies and the same types of 
configurations persist on both examined latices (6x6 and 8x8), confirming the 
stability of obtained results. 

Although we have presented here only the basic types of charge and magnetic su- 
perstructures (a more complete set will be given elsewhere) they clearly demonstrate 
an ability of the model to describe different types of charge and magnetic ordering. 
This opens an alternative route for understanding of formation an inhomogeneous 
charge/ magnetic order in strongly correlated electron systems. In comparison to pre- 
vious studies of this phenomenon based on the Hubbard and t — J model jH], the study 
within the generalized spin-one-half FKM has one essential advantage and namely 
that it can be performed in a controllable way (due to the condition [f^fia, H] = 0), 
and in addition it allows easily to incorporate and examine effects of various factors 
(e.g., an external magnetic field, nonlocal interactions, etc.) on formation of charge 
and magnetic superstructures. The work in this direction is currently in progress. 
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In summary, a combination of small-cluster exact-diagonalization calculations 
and a well-controlled approximative method was used to study the ground-state 
phase diagrams of the generalized spin-1/2 FKM with an anisotropic, spin-dependent 
on-site interaction between localized and itinerant electrons for Nf = L (0 < J < 1 
in D = 1 and < J < 4 in D = 2) and A^"/ < L (J = 0.5). For both cases it 
was observed that the anisotropic, spin-dependent interaction induces strong cou- 
pling between the localized and itinerant subsystems and that the magnetic phase 
diagrams of these subsystems coincide. In general, the spin-dependent interaction 
(for Nf = L) stabilizes the F and Fl state, while the stability region of AF phase is 
gradually reduced. The opposite effect on the F, Fl and AF phases has a reduction 
of Nf. For both Nf = L and Nf < L an extremely rich spectrum of charge and 
magnetic solutions has been found. In particular, we have observed various types 
of phase segregated, phase separated, striped, periodic and nonperiodic charge/spin 
distributions that clearly demonstrate strong cooperative effects of spin-dependent 
interaction on the ground states of the model. 

This work was supported by the Slovak Grant Agency VEGA under grant No. 
2/4060/04 and the Science and Technology Assistance Agency under Grant APVT- 
20-021602. Numerical results were obtained using computational resources of the 
Computing Centre of the Slovak Academy of Sciences. 
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Figure 1: The ground-state phase diagrams of the spin-one-half FKM extended by 
the Hubbard interaction between the itinerant electrons calculated for two different 
values of U on small finite clusters of L = 8, 10 and 12 sites. Below f/^^ the ground 
states are the ground-state configurations of the conventional spin-one-half FKM 
{Udd = 0). Above U^^ these ground states become unstable. The one-dimensional 
exact-diagonalization results. 
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O TTtTTTTttTtTTtTTi4.J,iiJ,iJ,4.iJ.ii4.iJ. 
T TTTTTTTTTTTi4.4.TTT4.J-iTTTJ-4.iTTT4.iJ, 

♦ TTTTTJ,itTTTTia.TTTTT4.iTTTTTJ.iTTiJ, 

• TTtTTTTtTTtTTTiiTTiiTTiitTJ-iTtiJ. 



Antiferromagnetic phases: 

■ (N^=6) tTTtTiJ,iJ,4.itTTtTJ.i4.J,iTTTttTJ.i4.4.i 

■ (N^=8) TTTt-LiiitTTTiiUTTTt-LiiitTTTiiU 

■ (N^=10) tTTiiiTTTTiiiTTTiiiTTTUiiTtTiU 

■ (N^=12) TTTJ-iTTTJ-iiTTJ-iiTTTJ-iTTTJ-iiTTiii 

■ (N^=14) tTTJ,iTTTJ,4.itTi4.itTTJ,iTTTJ,4.itTi4.i 

* (N^=14) tTiiiiiTTJ-iiiiTTiiTTTtniTTtTTU 

■ (N^=1 6) Tnj-TTiitTiJ.TTiitTiJ.TTiiTTiitTii 

■ (N^=18) tTiJ,TiJ,Tt4.itiitTJ.J,TtiTTiJ,tTJ.TT4.i 

■ (N^=20) TTiTTiTTJ-iTJ,iTi4.TTiTTiTTJ-iTJ,iTi4. 

■ (N^=22) TiitiiTiiTiiTiJ-TiTTiTTiTTJ-TtiTTJ. 

* (N^=24) TiitiTTitiiTiTTiTiitiTTitiiTitn 
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* (N^=26) TiTT4.TiTTiTJ,TiiTJ,TTJ-TiTiJ-TiTiiTi 
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■ (N =32) TiTJ,TiTiT4.TJ,T4.TiTiTJ,TiTiT4.TJ,TiT4. 



Figure 2: The /-electron ground-state phase diagram of the spin-one-half FKM 
extended by the spin-dependent interaction calculated for Nf = L {L = 32). The 
one-dimensional exact-diagonalization results. Inset: The d-electron ground-state 
phase diagram of the model calculated at the same values of model parameters. (■): 
the AF phase, (-I-): the F phase, (o): the FI phase. 
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Figure 3: The / electron skeleton phase diagram of the spin-one-half FKM extended 
by the spin-dependent interaction calculated for Nf < L, J = 0.5 and L = 24. (■): 
the AF phase, (-I-): the F phase, (o): the Fl phase. The one-dimensional exact- 
diagonalization results. 
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Figure 4: The / and d electron phase diagrams of the spin-one-half FKM extended by 
the spin-dependent interaction calculated for Nf < L, J = 0.5 and L = 48. (■): the 
AF phase, (-I-): the F phase, (o): the FI phase. The one-dimensional approximative 
results. 
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Figure 5: The / and d electron ground-state phase diagrams of the spin-one-half 
FKM extended by the spin-dependent interaction calculated for A'^^ = L [L = 36). 
The two-dimensional approximative results. 
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Figure 6: Typical examples of ground states of the spin-one-half FKM extended by 
the spin- dependent interaction obtained for A'^^ = L [L = 36). To visualise the spin 
distributions we use A for the up spin electrons and y for the down spin electrons. 
The two-dimensional approximative results. 



19 



(1) 



(2) 



(3) 



(4) 



(5) 



■ A A ■ 
A A A A 

■ AA • 



▼ ■ ■ 
Ay . 
■ At 

. . At 

... A 



■ ■ • ■ ▼ T 

.... AA 
▼▼ .... 
A A ... . 



A 
A 
A 
A 
A 
A 



(6) 



(7) 



(8) 



(9) 



(10) 



▼ ■ T 

. A . 



T 

A . A 



A . 
• ▼ 



A . 
■ T 



(11) 



(12) 



(13) 



(14) 



(15) 



A 
▼ 
T 

A 
▼ 

A 
A 
T 



. At ■ 

AAtt 
A AtT 

■ Ay ■ 



■ ▼ ■ 

A ■ A 



. A A 

A A 

A A A A AAA A A A 

T T 

T T ■ 

• • • • ▼▼ ■ ■ 

■ ■ ■ TT ■ ■ ■ 

■ ■ AA .... TTTTTTTT 



(16) 



(17) 



(18) 



(19) 



(20) 



■ T ■ T 

A . A . 



A . A . 
• T ■ T 



. . - . A A ▼ 

. . . AAy 

. . A A T 

. A A ▼ 

A A ▼ . 

A A y . . 

At . . . 

▼ ■ ■ ■ ■ 



A 
A A 



att 
att 
Att 
Att 
Att 
Att 
Att 
Att 



A . A . A . 
• T ■ T ■ T 

. . A . A . 



■ T ■ 

A . A 



T ■ ▼ 

. A . 



A . A . A . 
• T ■ T ■ T 



(21) 



(22) 



(23) 



(24) 



(25) 



AAA 
A T A 
AAA 
A T A 
AAA 
AtA 
AAA 

AtA 



. A 
T ■ 



A ■ A . A 
■ T ■ T ■ 



■ T ■ T 

A . A . A 



T ■ T • T ■ T 

. A . A ■ A - ^ 



A Att 
AAtT 

A Att 
A Att 
aatt 
A Att 
A Att 

AAtT 



A 
A 



A ■ A . A . A . 

■ T • T ■ T ■ T 

A . A . A . A . 

■ T ■ T ■ T ■ ▼ 

A ■ A . A . A . 

■ T ■ T ■ T ■ T 

A ■ A . A . A . 



(26) 



A 
T 



(27) 



(28) 



atta 
taat 
Att A 

T A At 
AttA 
T A At 

AttA 

T A At 



TTTTTTTT 

AAAAAAAA 

▼▼TTTTTT 

AAAAAAAA 



(29) 



T • A . T • A . 
TTTTTTTT 

. T • A . T ■ A 
AAAAAAAA 
T ■ A . T ■ A ■ 
TTTTTTTT 

. T ■ A . T ■ A 

AAAAAAAA 



(30) 



TTT ■ AAAt 
TTTAAA ■ T 
T . A A A A T T 
A A A A . TTT 
AAAtTTT ■ 
AA . tTTTA 
A A T T T ■ A A 
. T T T A A A A 



Figure 7: Typical examples of ground states of the spin-one-half FKM extended 
by the spin-dependent interaction obtained for Nf < L, J = 0.5 and L = 64. The 
two-dimensional approximative results. 
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